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terms of 1/z = 1/(toL ), which we test by a non-perturbative computation of quenched 
current matrix elements and energies, taking the continuum limit of lattice results. These 
are seen to approach the corresponding renormalization group invariant matrix elements of 
the static effective theory as the quark mass becomes large. We are able to obtain estimates 
of the size of the 1 /m-corrections to the static theory, which are also of practical relevance 
in our recent strategy to implement HQET non-perturbatively by matching to QCD in a 
finite volume. 
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1. Introduction 

It is expected that the dynamics of QCD simplifies in the limit of large masses of the charm 
and beauty quarks. For external states with a single heavy quark, transition amplitudes are 
expected to be described by an effective quantum field theory, the Heavy Quark Effective 
Theory (HQET) [1, 2, 3]. For details of the proper kinematics where this theory applies, 
and also for a guide to the original literature we refer the reader to reviews [4, 5]. To 
prepare for the following presentation, we only note that HQET applies to matrix elements 
between hadronic states, where these hadrons are both at rest and do not represent high 
excitations. HQET then provides an expansion of the QCD amplitudes in terms of l/m, 
the inverse of the heavy quark mass(es). The HQET Lagrangian of a heavy quark is given 
by 1 

£rqet(x) =Mx) \D + m-^n 2 -^a-B]Mx) + ■■■ > (1-1) 
L 2m 2m J 

where the ellipsis stands for higher-dimensional operators with coefficients of 0(l/m 2 ). 
Following power counting arguments, this effective theory is renormalizable at any finite 
order in l/m [6, 7]. A significant number of perturbative matching computations have been 



1 We write here the velocity-zero part, since non-vanishing velocities will not be relevant to our discussion. 
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carried out (see [8, 9] and references therein), in order to express the parameters of HQET 
(m, w S pi n , . . .) in terms of the QCD parameters. The very possibility of performing this 
matching reflects good evidence that HQET does represent an effective theory for QCD. 

Nowadays, HQET is a standard phenomenological tool to describe decays of heavy- 
light hadrons and their transitions in terms of a set of hadronic matrix elements, which 
are usually determined from experiments. Its phenomenological success is illustrated, for 
example, by the determination of the Cabibbo-Kobayashi-Maskawa matrix element V c ^,: 
its value determined from inclusive b — > c transitions agrees with the one extracted from 
exclusive ones [10, 11] and HQET enters in both determinations. Similarly, HQET hadronic 
matrix elements, such as Ai oc oJk in (B| — Y> h D 2 ^h|B), extracted from different experiments 
tend to agree [11, 12, 13]. As a small caveat concerning these phenomenological tests, we 
note that some of them involve both the beauty and the charm quark, and one may suspect 
that the effective theory is not very accurate for the latter. 

Additional, independent tests of HQET are thus of both theoretical and phenomeno- 
logical interest. In [14] some Euclidean correlation functions, which are gauge invariant and 
infrared-finite, were studied in perturbation theory. There it was verified at one-loop order 
that their (large-distance) m — > oo asymptotics is described by the correlation functions of 
the properly renormalized effective theory. This comparison was performed after separately 
taking the continuum limit of the lattice regularized effective theory and of QCD. 

In general, lattice gauge theory calculations allow a variation of the heavy quark mass 
and the performance of non-perturbative tests. However, if one is interested in the com- 
parison of QCD and HQET in the continuum limit, one has to first respect the condition 

m<-, (1.2) 
a 

and then to do an extrapolation to zero lattice spacing, a = 0. Given the present restrictions 
in the numerical simulations of lattice field theories, a direct comparison at large mass can 
only be done in a finite volume of linear size significantly smaller than 1 fm. 

Before discussing this further, we note that in the charm quark mass region the con- 
tinuum limit can also be taken in a large volume [15, 16, 17]. In particular, a recent 
study concentrated on the decay constant Fd b , where many previous estimates found evi- 
dence for large 1/m-corrections. After taking the continuum limit and non-perturbatively 
renormalizing the lowest-order HQET, ref. [18] finds only a rather small difference between 
lowest-order HQET and the QCD results. Although this investigation was restricted to 
the quenched approximation, it provides some further evidence of the usefulness of HQET 
- maybe even for charm quarks. For related work, see [19, 20, 21, 22, 23, 24, 25, 26, 27, 
28, 29, 30, 31] and references therein. 

In this paper we study the large-mass behaviour of correlation functions in a finite 
volume of size L x L 3 , with Dirichlet boundary conditions in time and periodic boundary 
conditions in space, i.e. we work in the framework of the QCD Schrodinger functional 
[32, 33]. Keeping all distances in the correlation functions of order L, the energy scale 1/L 
takes over the role usually played by small (residual) momenta, and at fixed L, HQET can 
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be considered to be an expansion of QCD in terms of the dimensionless variable 



1 



1 



(1.3) 



z 



ML ' 



where we take M to be the renormalization group invariant (RGI) mass of the heavy 
quark (see eq. (3. 6)). 2 With the choice L ~ 0.2 fm, today's lattice techniques allow us to 
increase z beyond z = 10, while the continuum limit can still be controlled well [34, 35]. 
One is thus able to verify that the large- z behaviour is described by the effective theory, 
which is the primary purpose of this paper. Of course, the coefficients a« in expansions 
<E> = ao + a\jz + a,2 /z 2 + . . . are functions of AL (with A the intrinsic QCD scale) such that 
aj(AL) — > ci x (AL) 1 and a^KVj/z 1 — > q x (A/M) 1 as AL — > oo. Since we only work at 
one value of L, the dependence on AL will be suppressed in the rest of this paper. 

After a discussion of the QCD observables under investigation (section 2), we give their 
HQET expansion (section 3) and confront them with Monte Carlo results at finite values 
of 1/z (section 4). In our conclusions we also discuss the usefulness of our results for a 
non-perturbative matching of HQET to QCD [36]. 

2. Observables 

We introduce our observables starting from correlation functions defined in the continuum 
Schrodinger functional (SF) [32, 33]. For the reader who is unfamiliar with this setting, we 
give a representation of these correlation functions in terms of operator matrix elements 
below. 

We take a T x L 3 geometry with T/L = 0(1) fixed. The boundary conditions are 
periodic in space, where for the fermion fields, and only for those, a phase is introduced: 



In the numerical investigation of section 4, we will set T = L and 6 = 0.5. Dirichlet 
conditions are imposed at the xq = and xq = T boundaries. Their form as well as the 
action are by now standard [32, 33, 37] and we do not repeat them here. Multiplicatively 
renormalizable, gauge-invariant correlation functions can be formed from local composite 
fields in the interior of the manifold and from boundary quark fields. Boundary fields, 
located at the xq = surface are denoted by £, £ and create fermions and anti-fermions. 
Their partners at xq = T are (',('■ We shall need flavour labels, "1" denoting a massless 
flavour, "b" a heavy but finite-mass flavour and "h" the corresponding field in the effective 
theory. The heavy-light axial vector and vector currents then read 



2 We continue to use m as a generic symbol for the quark mass when its precise definition is irrelevant, 
sometimes even not distinguishing the renormalized mass from the bare one. However, when precisely 
defined functions of the quark mass are considered, only the properly defined M is used. 



ip(x + kL) = e ie V(x) , xj}{x + kL) = ^(x)e~ ie , k = 1,2,3. 



(2.1) 



(2.2) 
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2.1 Correlation functions 

In our tests we shall consider the correlation functions 

/a(* ) = -\\ d 3 yd 3 z <A)(*)Cb(y)75Ci(z)} , 



kv(x ) = -\Y.j d3 y d3z W*Kb(y)7*Ci(z)) , 



as well as the boundary-to-boundary correlations 
1 



fi = - 



2L 6 



d 3 yd 3 zd 3 y ' dV <C;(y')75C b (z') C b (y)7 5 Ci(z)> 



= -«P E / d 3 yd 3 zd 3 y'd 3 z' <c((y / )7 fc C b (z / ) Cb(y)7.Cl(z)> 

They are illustrated in figure 1. 
xq = T 

A (x ) V k (x ) 



(2.3) 
(2.4) 

(2.5) 
(2.6) 



Cb 75 cl 






x = 



Figure 1: The Schrodinger functional correlation functions /a, &v and f\. For fci, in the rightmost 
diagram 75 is replaced by j/.- 

Let us discuss the function /a(xo) in some detail. It describes the creation of a (finite- 
volume) p = heavy-light pseudoscalar meson state, \(pb(L)), through quark and antiquark 
boundary fields which are separately projected onto momentum p = 0. "After" propagating 
for a Euclidean time interval xq, the axial current operator Ao initiates a transition to a 
state, |f2'(L)), with the quantum numbers of the vacuum. The correlation function /a(xo)> 
taken in the middle of the manifold, can thus be written in terms of Hilbert space matrix 
elements, 

/ A (T/2) = Z-\Sl(L)\ A \B(L)) = Z-\Sl(L)\ti{L)) , \B{L)) = e~™/ 2 \ML)) , (2.7) 



where EI is the QCD Hamiltonian and 

z = (n(L)\n(L)) , \Q(L)) = e-™i 2 \ m {L)) . 



(2.8) 



Here, \ipo{L)) denotes the Schrodinger functional intrinsic boundary state. It has the quan- 
tum numbers of the vacuum. Because the definition of our correlation functions contains 
integrations over all spatial coordinates, all states appearing in our analysis are eigenstates 
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of spatial momentum with eigenvalue zero. The operator e _TH//2 suppresses high-energy 
states. When expanded in terms of eigenstates of the Hamiltonian, \Q(L)) and \B(L)) are 
thus dominated by contributions with energies of at most AE = 0(1/ L) above the ground 
state energy in the respective channel (recall that we take T = 0(L)). This explains 
why, at large time separation xq > 0(l/m), HQET is expected to describe the large-mass 
behaviour of the correlation function, also in the somewhat unfamiliar framework of the 
Schrodinger functional. 3 

Equations similar to the above hold for ky, one only needs to replace pseudoscalar 
states by vector ones. Finally, the boundary-to-boundary correlator may be represented as 

f 1 =Z- 1 (B(L)\B{L)). (2.9) 

Since the boundary quark fields are multiplicatively renormalizable [38], this holds 

also for the states \ipo(L)) and \ipb(L)}. 
It now follows that the ratios 

YML,M) = f -^, y v( L,M).-M|a, R{LtM) ^-^m (2.10) 



are finite quantities when we adopt the convention that A^,V^ denote the renormalized 
currents. As is immediately clear from the foregoing discussion, 

PS( ' } = || ML)) B(L)) || (2 - U) 

(or Yy(L, M)) becomes proportional to the pseudoscalar (or vector) heavy-light decay 
constant as L — ► oo. We shall study the large-M behaviour of these quantities at fixed L 
in the following sections. 

For the same purpose we define effective energies 



r PS (L,M) = -A-ln[f A (x )} 

T V (L,M) = -Ai n [& v ( Xo )- 
dxo 



x =T/2 



x =T/2 



/a(T/2) ' [2A2) 



ky(T/2) 



(2.13) 



These may be written in terms of matrix elements of the Hamiltonian. E.g., we have 

(B\L)\a\B(L)) (ii(L)|H|ii'(L)> 

r PS (L,M)= wmm) - {(1(L)mL)) ■ 



with \B'(L)) = A t |^(L)) and (B'(L)\B(L)) = {Sl(L)\SV \L)) . Expanding in terms of 
energy eigenfunctions, one sees immediately that r PS (L,M) = YliPi E< B ~ J2i a i E n\ 
where Eg are the (finite-volume) energy levels in the heavy-light pseudoscalar meson 

(i) 

sector and E^ those with vacuum quantum numbers. The coefficients (5i,ai have a strong 
dependence on i, which labels the energy excitations; states with E^ — E^ » 1/L are 



3 More generally, HQET will apply to correlation functions at large Euclidean separations. 
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suppressed exponentially in the sum. For z = ML S> 1, the effective energy Tps(L,M) is 
hence expected to be given by HQET. 

To summarize, Yps, Yy and R are (ratios of) matrix elements between low-energy 
heavy- light and vacuum-like states. HQET at order (l/m) n should describe them up to 
corrections of the order of \/{ML) n+1 . The energies T, eqs. (2.12) and (2.13), share the 
same property. It is then possible to test HQET by studying the large- z asymptotics of 
these observables! 



3. Large-z asymptotics: effective theory predictions 

We now turn our attention to the effective theory predictions for the observables introduced 
above. 



3.1 Current matrix elements 

At the classical level it is expected that they can be described by a power series in 1/z with 
z = ML. This has been checked explicitly in [14]: the expansion in 1/z is asymptotic, 
non-analytic terms are of the type e~ z and are thus very small for, say, z > 4. The leading 
term in the expansion for each of the correlation functions in eqs. (2.3)— (2.6) is given by 
exactly the same expressions, evaluated with the simple replacement ifa, — > Vh etc- an d by 
dropping all terms of 0(l/m) in the action associated with eq. (1.1). This corresponds 
to the static limit, where the heavy quark does not propagate in space. We denote the 
corresponding observables (in the effective theory) with a superscript "stat" , e.g. 

A0roWX tat 0£o), 

and also introduce 

XiL) = ffWl (SJ) 

that is easily seen to be mass independent from the form of the static propagator (both in 
the formal continuum theory and in lattice regularization at finite lattice spacing [37]). An 
example for the correspondence of the effective theory and QCD at the classical level is 

X(L) = lim Y PS (L,M) = lim Y V (L,M). (3.2) 

z — >oo z — >oo 

In the quantum theory there are logarithmic modifications of such relations. Of course, 
they have their origin in the scale dependent renormalization of the effective theory. For 
example, due to the renormalization of the axial current in the effective theory, the renor- 
malized ratio 4 



X R (L, f i) = Z s ^ t (fi)X haTC (L) (3.3) 

4 One should not identify fx = 1/L here, as it was done in the computation of the scale dependence of 
the static axial current in [39]. Note also that in the ratio X other renormalization factors cancel in the 
effective theory, just as they do in QCD. 
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depends logarithmically on the chosen renormalization scale \x. It further depends on the 
chosen renormalization scheme, but the renormalization group invariant 

X RG1 (L) = lim { [ 2b g 2 (n) } " 7o/(26o) X R (L, M )| = Z RGl X harc (L) , (3.4) 

6 ° = 7^)2> 70 = -^) (JVf = 0), (3.5) 

does not. The renormalization constant Z RG \ is computable in lattice QCD [39]. Above, 
bo and 70 are given for the case of a vanishing number of flavours as is appropriate for the 
quenched approximation, which we will employ in the following section. 

The large- z behaviour of the QCD observables is then given by the corresponding 
renormalization group invariants of the effective theory, together with logarithmically mass 
dependent functions that will be called C below. As arguments of these functions we choose 
the renormalization group invariant mass M (in units of the A parameter), since this can be 
fixed in the lattice computations without perturbative uncertainties: the relation between 
the bare quark masses in the lattice regularization, which we use [40], and M has been 
non-perturbatively computed in [41, 42, 35]. The scheme independent M describes the 
limiting behaviour of any running mass m(/x) for large /i via 

lirn{[26o5 2 (^]" d ° /(2M ^)}=M, d = (N { = 0) . (3.6) 

The following predictions are then obtained (see also section 5.1 of ref. [39]): 

Y PS (L, M) M ~°° C PS (M/A m ) Ar GI (L) (l + 0(1/*)) , z = ML , (3.7) 
Y W (L, M) M ~°° CV (M/A m ) X RGI (L) (l + 0(1/*)) , (3.8) 
R(L,M) M ~°° C PS /v(M/A M g) (l + 0(l/z)). (3.9) 

Here, the function Cps(M/Aj^g) has the asymptotics 

-(^)-(,4)-^{ 1+ o(^^)} ; «,„, 

to within this accuracy, Cy shares this asymptotic behaviour. In practice, the functions 
Cx(M/Aj^g), X = PS,V, PS/V, are obtained by solving the perturbative renormalization 
group equations along the lines of section 5.1 of ref. [39], where more details can be found. 
In particular we always use the four-loop perturbative approximation of the /^-function [43] 
and the three- loop approximation to the anomalous dimension 7 of the currents, which has 
recently been computed [9]. Taking the n-loop approximation to 7, the resulting relative 
error in the functions Cx is of order a n with a evaluated at a scale of the order of the 
heavy quark mass. Explicit expressions for the functions Cx are given in appendix B. 

In order to study whether their a n -error is a limiting factor, we first note that changing 
the order from three-loop to four-loop in the /3-function makes only a tiny difference. Next 
we plot Cx, X = PS, PS/V, using 7 at n loops in figure 2. We observe a reasonable 



-7- 



0.2 



0.4 



0.6 1/z 



t 1 r - 



1 -loop 7 - 

2-loop 7 

3-loop 7 




J PS/V 

1.1 



1.05 - 



1 - 



h 1 1 1 1 r 

1-loop 7 

2-loop 7 

3-loop 7 




0.05 



0.1 



Figure 2: The functions Cps and Cps/v evaluated in different approximations of perturbation 
theory. In relating A^jg/Af to z — ML, we used [41] A^g-r = 0.602 and L/r = L max /(2r ) = 
0.359. (Setting tq = 0.5 fm, this actually corresponds to L m 0.2 fm.) 



behaviour of the asymptotic perturbative series and take half of the difference between the 
two-loop and three-loop approximations for 7 as our uncertainty. This uncertainty will be 
almost negligible with respect to our statistical errors. Note that without the three-loop 
computations of [9] such a statement would not have been possible. 

For completeness we remind the reader that eqs. (3.7)-(3.9) are afflicted by the usual 
problem of identifying power corrections. Asymptotically, at large M, the higher-order 
logarithmic (perturbative) corrections in Cx dominate over the l/z n power corrections. 
However, as just discussed, in the interesting range of z, the perturbative corrections are 
under reasonable control and it makes sense to investigate the power corrections if they 
are significantly larger. Note that this problem is not present if the effective theory is 
renormalized non-perturbatively [36]. 

3.2 Energies 

Next we concentrate on the energies T. Clearly, because of the term mt/; h (x) , ip\ i (x) in the 
Lagrangian (1.1), they grow (roughly linearly) with the mass. For the case of hadron 
masses, mass formulae have been written down [44]; for instance, the one for the B-meson 
mass reads: 

m B = m + A+^-(X 1 +d B X2) + O (l/m 2 ) , d B = 3 (3.11) 
2m 

(and the same formula holds for the B*-meson except that d B — ► d B * = — -0- The matrix 
element Ai = ^Wkin(B| — ip^D 2 ^^) was mentioned already in the introduction, and 
^ 2 = 23i" Ws P in ^l — ^h 17 " B^h|f3). Some cautioning remark concerning the above formula 
is in order. It suggests that the binding energy A = m B — m + 0(l/m) may be obtained as 
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a prediction of HQET. However, there is no unique non-perturbative definition of the mass 
m, which should be subtracted. As a consequence, A has an ambiguity of order Aqcd, 
which may then also propagate into a significant ambiguity in Ai extracted from eq. (3.11). 

We nevertheless start our discussion of the effective energies T from a trivial general- 
ization of eq. (3.11): 



r (L,M) = 1 r^(L.M) + :\\: v (L. U) 



= m + A(L) + — Ai(L) + 0(l/m 2 ) 



(3.12) 



where Ai(L) again summarizes the effect of the V> h D 2 ?/;h-perturbation to the static action. 
We have cancelled a A2-term by considering the spin-averaged combination of the energy. 
While one ought to be careful with the interpretation of sub-leading terms in 1/M, the 
large-mass behaviour of T av (L, M) is given by 



Lr av (L,M) 



M^oo 



C mass {M/A m ) xz + O ((l/z)°) , 



with (mj^g(m*) = m* 



C mass (M/A^) = 



M 



mg 



(3.13) 



(3.14) 



and niQ being the pole mass. Here the first factor on the right-hand side is known very 
precisely in perturbation theory (up to four loops [45, 46]), but it is well known that the 
perturbative series for the second factor is not very well behaved and even the three-loop 
term [47] is still rather significant. We will discuss this uncertainty in C mass together with 
the numerical results in the following section. 5 

Let us now consider the sub-leading terms in eq. (3.12). Taking energy differences, m 
drops out and e.g. A(L) — A(L') can be computed unambiguously from the static effective 
theory. It does not depend on any convention adopted for m in eq. (3.12). In our numerical 
computations, however, we investigated only one value of L. As an example of another 
observable unaffected by the ambiguity in m, we therefore look at the combination 



1{L,M) = \ 



■ f' A (T/4) f A (T/2) k' y {T/A) k'yjT/2) 
,/a(T/4) /a(T/2) + fcv(T/4) fcv(T/2) 



The HQET prediction for this quantity is 



E(L, M) = S stat (L) + j- H kin (L) + O (l/z 2 ) , 



with 



^stat(^) — L 



(/i taty (T/4) 

/rtat (T/4) 



(/i taty (r/2) 

/rtat (T/2) 



(3.15) 



(3.16) 



(3.17) 



5 We note that eqs. (3.13) and (3.14) are the only ones where rriQ enters the coefficient functions relating 
the RGI matrix elements of HQET to the QCD observables. In all other cases, rriQ has been eliminated 
and only M appears. Thus, there is no particular reason to expect that any of the perturbative expressions 
(for the various anomalous dimension functions) is badly behaved. 
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denned in the static effective theory. Since it is an energy, S does not require any renor- 
malization. The first-order correction in 1/z is given entirely by matrix elements of 
ip h (x)D 2 ipi 1 (x)- Its coefficient is fixed using the reparametrization invariance of the ef- 
fective theory, first discussed in ref. [48]. To see this, one considers the matrix elements 
of ip h (x)Y) 2 ip\ 1 (x) to be computed in dimensional regularization. Then reparametrization 
invariance is valid, and the coefficient of the operator renormalized by minimal subtraction 
is the inverse MS quark mass [49, 50], whose renormalization scale dependence cancels 
against the one of the matrix element. From eq. (3.6) it hence follows that the prefactor of 
the renormalization group invariant matrix element of ip h (x)Y) 2 ip\ 1 (x) is 1/(2M), as used in 
eq. (3.16). In that equation, E^ in (L) is the matrix element made dimensionless by a factor 
L. In the following, the exact expression for £khi(L) will be irrelevant. It rather suffices to 
know that it does not depend on the mass. 

We point out that, in writing down the quantum mechanical representation of S, 
states play a role for which only the operator i s effective to suppress higher energy 

contributions, while in our other observables e -™/ 2 suppresses high energies. For the 
quantity H, HQET is thus expected to be accurate only at larger values of z. This variable 
should roughly be a factor of 2 larger than for the other observables. 

Finally, in the difference 

A r (L, M) = r PS (L, M) - F V (L, M) (3.18) 

the lowest-order term that contributes in the effective theory is iph( x )< J ' B^h(x). With 
C sp in constructed from the anomalous dimension of this operator in the effective theory, 
7 spin [51, 52], we thus have 

-v^spin / j \ 

LA r (L, M) M ~°° C spin (M/A m ) -^J^ (l + 0(l/z)) , (3.19) 
where the leading asymptotics of the function C sp i n (M/Aj^g) is of the form 

CWM/An,) ~ (M—j | 1 + ( ln( M/A w ) )} (320) 

with the universal coefficient (cf. appendix B) 

(3 ' 21) 

and -^rgi being again a renormalization group invariant matrix element, defined in the 
static effective theory. 6 The relative perturbative uncertainty of C sp i n is 0(q 2 ), since here 
the three-loop anomalous dimension is not known. As illustrated in figure 3, the difference 
between the one-loop anomalous dimension and the two-loop one is tiny. Since this may, 
however, be accidental rather than representing the behaviour of the series, we shall take 
the size of the three- loop term met in Cps, figure 2, as our uncertainty. A better estimate 
of this uncertainty would require the knowledge of 72 Pm - 

6 Alternatively, one can also construct a difference of squared energies as the product L 2 r av Ar, which 
then behaves as c7 mag X^Qj (1 + 0(1/2)) for M — > 00, where C mag is introduced in appendix B as well. 
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Figure 3: The function C sp i n evaluated in perturbation theory. 



4. Quantitative tests of the effective theory 

We have tested the predictions of the effective theory applied to the quenched approximation 
of QCD by evaluating the observables for one geometry, namely T = L = Lq, and for 
9 = 0.5. The use of the quenched approximation should not be worrying in this context, 
since it is of course used both in the effective theory and in QCD. Furthermore, although 
we set the light quark mass to zero, 1/Lq provides an infrared cutoff and there are no 
singularities in the chiral limit. 

The numerical simulations are done on lattices with various resolutions a/Lo followed 
by a continuum extrapolation. In all cases, Lq is fixed by imposing 



where g(L) is the renormalized coupling at length scale L in the SF scheme [53]. It is known 
that Lq ~ 0.2 fm [41, 36]. The (purely technical) reasons for the precise definition (4.1) are 
detailed in ref. [36]. Table 1 of ref. [35] lists the bare coupling go for each resolution Lo/a, 
and this reference also explains how the bare quark masses are fixed to ensure a massless 
light quark and a prescribed value M for the heavy quark. 

4.1 Results in the static approximation 

Some of the leading-order terms in the HQET expansions described in the previous section 
are known exactly due to the spin symmetry [54, 55], but for two of the expansions we have 
computed the non-trivial leading order from a lattice simulation in static approximation. 

The first one is the matrix element of the time component of the axial current, X R qi(L) 
in eq. (3.4). Indicating explicitly also the dependence on the bare coupling go and the lattice 
spacing a, it is given by 

X RGl (L )= lim Z RG1 {g )X{g ,L /a). (4.2) 

a/Lo^0 

From the relation of X to the very definition of the renormalization factor Zrgi , as detailed 
in ref. [39], one obtains the explicit form 



g 2 (L /2) = 1.8811 



(4.1) 




x [X] 



tree— level 



(4.3) 
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/V RGI 

-1.278 
-1.28 
-1.282 
-1.284 

""stat 

0.12 
0.1 
0.08 
0.06 

0.002 0.004 0.006 (a/L) 2 

Figure 4: Continuum limit extrapolations of Xrgi(^o) (top) and 2 sta t(io) (bottom). The fits 
adopt simulation data generated with the HYP-link static quark action, see text. The a/L inde- 
pendent part of the error of Xrgi is not included in the graph. 

for the right-hand side of eq. (4.2). All quantities that enter the above expression have been 
defined in the latter reference. For its numerical evaluation we take the non-perturbatively 
O(a) improved action of [40] for the light quarks and the improved discretizations of [56] 
for the heavy quark. We note in passing that, using the data of [39], we first evaluated this 
quantity with the Eichten-Hill action for the heavy quark [1]. However, this resulted in an 
order of magnitude larger error for Xrgi at Lq/cl = 32. 

The limit a/Lo — » is taken by a linear fit in (a/Lo) 2 as illustrated in the upper 
diagram of figure 4, referring to the data set from a simulation with the static quark action 
built from HYP-links [57, 56]. We quote the result from a fit with Lq/cl > 16 as our 
continuum result, 

X RG i(L ) = -1.281(9), (4.4) 

which also receives an error contribution from the (regularization independent) factor 
$RGl/$(l/-£o) in eq. (4.3) [39]. Including all points with Lo/a > 12 in the fit yields a 
compatible continuum value with smaller error. In the same way we obtain (cf. eq. (3.17) 
and the lower diagram of figure 4) : 

Hstat(A)) = 0.11(1) . (4.5) 
4.2 Results at finite M and comparison 

The finite-mass (quenched) QCD observables are obtained from similar extrapolations of 
lattice results at finite a/L. However, as the variable z is increased, the quark mass in 
lattice units grows at a given resolution a/L. A perturbative computation [14] as well 
as our non-perturbative study suggest that O(a) improvement may be trusted only below 
a certain value of the quark mass in lattice units. It is therefore necessary to impose a 
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-1.6 



0.002 0.004 0.006 (a/L) 2 

Figure 5: Continuum extrapolations for some ^-values spanning the entire range of z. The linear fit 
in (a/L) 2 is shown by the full line and extended also to those values of a/L that did not participate 
in the fit by a dotted line. 

cut on the quark mass. For a given z, this cut translates into estimates of the coarsest 
lattice resolutions beyond which the lattice data are to be omitted from the continuum 
extrapolations. As in ref. [35] we impose aM < 0.6 and follow this reference on all other 
details of the lattice computations as well as the extraction of the continuum limits. For 
illustration we just show the continuum extrapolations of Yps(L,M), L = Lq, at selected 
values of z in figure 5. In addition we mention a few features equally true for the continuum 
limit extrapolations of the other observables, which enter our investigation but are not 
shown in figures. 

• The slopes in a/L are rather small. 

• The error in the continuum limit grows with z, because less lattices can be used in 
the extrapolation at large z. 

• The continuum limit is compatible with the values at the smallest two lattice spacings. 
Its error is conservative. 

In appendix A we collect our numerical results both at finite lattice spacing and in the 
continuum limit. 

We are now ready to compare the continuum results with the predictions of HQET. 
Before that, we remind the reader that energies of order 2/L still contribute significantly 
to our observables. It is thus possible that the 1/M-expansion breaks down earlier in our 
finite- volume situation than it does in large volume. However, our numerical results do not 
give any indication of such a behaviour. 

Let us start with the current matrix elements. The prediction for the matrix element 
of the axial current is Yps(L,M)/Cps = Xpg,i(L) + 0(l/z). In this combination, plotted in 
figure 6, the perturbatively computed coefficient Cps(M/Aj^g) compensates a significant 



Linear Quadratic 
Quantity ao a\ ao ai a<i 



z-range: 3.0-13.5 



yps/Cps 


-1.281(7) 


0.64(8) 


-1.0(2) 


IV/Cv 


-1.281(7) 


-0.63(9) 


0.3(2) 


-R/ Cps/v 


1.0 


-0.89(1) 


1.06(3) 


J-^sv / (.Z Cmass) 


1.0 


0.42(3) 


0.14(10) 


LAr/C S pi n 


0.0 


-1.69(6) 


0.8(2) 




0.109(8) 


0.7(1) 


-0.8(3) 



z-range: 5.15-13.5 



yps/Cps 


-1.277(7) 


0.45(4) 


-1.281(8) 


0.7(1) 


-1.3(5) 


Yy/Cy 


-1.281(9) 


-0.59(6) 


-1.281(7) 


-0.6(1) 


0.2(6) 


R l Cps/v 


1.0 


-0.722(7) 


1.0 


-0.91(1) 


1.18(5) 


-^r av / (z C m ass) 


1.0 


0.44(2) 


1.0 


0.41(6) 


0.2(3) 


LAr/C S pi n 


0.0 


-1.56(4) 


0.0 


-1.62(6) 


0.4(2) 




0.112(8) 


0.54(7) 


0.110(9) 


0.6(2) 


-0.5(7) 



Table 1: Fit parameters describing the z-dependences of the form a + a\/z + a 2 /z 2 . Where the 
constant term, ao, is listed without an error, it is constrained to the prediction of the static effective 
theory. (In case of the entire z-range, only the quadratic fit results are given.) All these fits have 
an acceptable goodness-of-fit. 



part of the mass dependence of Ypg. The finite-mass Ips/Cps is obviously well compatible 
with approaching the static result, eq. (4.4), as 1/z ^ 0. To quantify the deviations 
from the static limit at finite z, we fit all data points, including 1/z = 0, to first- and 
second-order polynomials in 1/z: 

Yps . a i u «n 

— = a + — + . . . . (4.6) 
Gps z 

In these fits also the uncertainty in Cps is taken into account (i.e. half of the difference 
between Cps evaluated with the two- loop and three- loop anomalous dimension). We per- 
form them separately to the data in the range 1/z < 0.2, which means masses around 
the b-quark mass and higher [36], and over the whole range, see table 1. Comparing a\ 
obtained from the quadratic fit over the whole range (z = 3.0-13.5) and the linear fit 
for 1/z < 0.2 (z > 5.15), the change is not so small. This indicates that a precise iden- 
tification of the first-order correction is not possible with our precision and range of z. 
However, the rough magnitude of a\ ~ 0.6 can be inferred and, more importantly, it is 
clear that the overall magnitude of 1/z-corrections is reasonably small. It is also relevant 
to remember that eq. (4.6) is only an approximate parametrization of the z-dependence, 
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Figure 6: The ratio >ps(L,M)/Cps(M/Aj^g) with Cps using the two- and the three-loop anoma- 
lous dimension of the static axial current. Error bars do not contain the perturbative uncertainty 
in Cps • The fits shown refer to the three- loop evaluation of Cps and include the result for Xrgi in 
the static limit. (Points using only the two-loop anomalous dimension are slightly displaced here.) 
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Figure 7: As in figure 6. The fits use the evaluation of Cps/v with the two-loop matching coefficient 
between HQET and QCD; they are constrained to the prediction of the static effective theory. (On 
the level of perturbative orders, two-loop matching belongs to the three-loop anomalous dimension 
of the currents, cf. appendix B). 



since the renormalization of the higher-dimensional operators in the effective theory will 
introduce logarithmic modifications of the simple power series. It is thus conceivable that 
these logarithms account for some of the curvature seen in figure 6. 

For the ratio of matrix elements R(L,M), eq. (3.9), the lowest-order term in the 1/z- 
expansion is fixed to be 1 by the spin symmetry of the static theory. As is reflected by 
figure 7, the results at finite 1/z are well compatible with this, if the function Cpg/y is 
evaluated including at least the two-loop anomalous dimensions. Fits to R/Cpg/y are 
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Figure 8: Quadratic fit over the entire z-range of the combination LT av /(zC masa ), which is con- 
strained to approach 1 in the static limit. It employs the three- loop relation between the pole and 
the MS mass; the results using this relation to only two loops are also depicted for comparison. 
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Figure 9: Fits in 1/z constrained to the static result and referring to the two- loop evaluation of 



Csi 



The linear fit is only based on the heavier quark mass points with z = 5.15-13.5, while the 



quadratic one (i.e. also allowing for a l/z 2 -term) includes all points. 



performed in complete analogy to eq. (4.6). The corresponding parameters of table 1 are 
again of order unity. In that table we also include the parameters of the analogous fits to 
the quantity Yy/Cy. 

Turning our attention to the effective energies, introduced in section 2.1, a first rough 
test of HQET is the behaviour of the spin-averaged energy T av . Figure 8 confirms the 
expectation that the combination LT av /(z C mass ) approaches 1 up to 1/z-corrections as 
1/z — ► (see eq. (3.13)). Note that the (l/z) 2 -corrections in T av are very small, which is of 
particular interest for the computation of the b-quark mass in static approximation [36, 35]. 
There the quantity r av was used in order to non-perturbatively match the quark mass of 
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the effective theory to the one in QCD. The dominant error in the final estimate for the 
quark mass, Mb, is expected to originate from this matching and is of order My, x (1/z) 2 . 
From the values of 02 and z\, = m^L ~ 5, we may estimate the relative error to be roughly 
of the order of 0.2 x (1/5) 2 ~ 1%. As seen from the figure, this conclusion is not very 
much affected either by the perturbative uncertainty visible in the graph as a difference 
between the two- loop and three-loop approximations for the pole mass tuq in eq. (3.14); 
the l/z 2 -curvature is not very different. 

The spin splitting Ap(L, M), which vanishes in the static limit, is displayed in fig- 
ure 9. It is in good agreement with the HQET prediction but exhibits a rather large 
1/z-coefficient. 7 

Finally, we successfully test eq. (3.16) in figure 10. Recall that a simple kinematical 
consideration leads one to expect the l/z-expansion to only be applicable at smaller values 
of 1/z for this observable (cf. section 3.2). On the other hand, reparametrization invariance 
allows to exclude logarithmic modifications of the 1/z-term; in contrast to our other tests, 
figure 10 does not involve any perturbative conversion factor. 
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Figure 10: Linear and quadratic fits of the observable S as in figures 6-9, where the static data 
point, S s tat, is included in the fits. 

5. Conclusions 

All of the comparisons of QCD observables with the predictions of HQET discussed in this 
work represent tests of the effective theory, which it passed successfully. As a significant 
improvement of earlier non-perturbative (lattice) tests, these comparisons are performed 
after first taking the continuum limit of the non-perturbatively renormalized quantities. 

7 We note that the slope X^^(L) is computable in the effective theory in a very similar way to Xrgi, 
because the associated operator ip h a ■ B iph does not mix with any other; its renormalization may be 
computed non-perturbatively using the methods of ref. [39]. The comparison of the result to the data at 
finite mass would presumably be limited in precision by the present perturbative uncertainty in C sp m- This 
limitation is likely to apply to the (large- volume) mass splitting between the B*- and the B-meson as well. 



-17- 



Thus, worries that 0(a) effects, in particular those that grow with the quark mass, may 
afflict the large-mass behaviour in QCD are removed. 

For a precise judgement of figures 6-10 it is important to be aware of the level of 
precision of renormalization and matching. The renormalization problem of the static axial 
current has been solved non-perturbatively in [39] and with this information all points at 
1/z = are known without any residual perturbative uncertainty. For the quantity E 
shown in figure 10, this is also true at finite 1/z. However, in general we need to know the 
functions Cx, which relate the QCD observables at finite mass M (and thus finite z = ML) 
to the renormalization group invariants of the effective theory, such as -Xrgi m figure 6. 
While the latter are unambiguously defined and have been computed non-perturbatively, 
the former are known only in perturbation theory and have errors of order a n , with a 
evaluated at the scale of the heavy quark mass. We have discussed these errors. On a 
phenomenological level they are under control due to the computations [58, 9], except for 
Cspin, where a next-to-next-to-leading-order computation is not yet available. Still, one 
should remember that - strictly speaking - the isolation of 1/m-corrections by looking at 
the difference to the leading-order HQET result is only possible when Cx is known non- 
perturbatively: parametrically, 0{[a{m)] n ) corrections are always larger than 0(l/m). 
Nevertheless, our results in figures 6-10 are compatible with the z~ n power corrections 
dominating over the perturbative ones in the accessible range of z. Fitting them by a 
simple polynomial in 1/z, the coefficients turn out to be of order 1 as naively expected. 
Therefore, at the b-quark mass, which corresponds to z^ = m^L « 5 for our value of 
L = Lo ~ 0.2 fm, the effective theory is very useful. 

Beyond the general interest of providing a non-perturbative test for HQET, these 
results are important for the programme of ref. [36]. There it was suggested to determine 
the 1/mb-corrections to B-physics matrix elements from a simulation of HQET on the 
lattice. The coefficients Ck of the various terms in the HQET Lagrangian and of the 
effective composite fields are to be determined by matching HQET and QCD in a finite 
volume, and it was proposed to employ a value of L similar to Lq in this step. This value 
has to be large enough such that (i) z^ is in the range where HQET is applicable with 
small corrections, yet (ii) Lq has to be small enough to allow for the computation of the 
QCD observables at m = with small a-effects. From table 1 we conclude that L = Lq is 
indeed promising. In fact, in the previous section we roughly estimated that the correction 
to the static limit computation of the RGI b-quark mass [34, 36] is only of the order of 1%. 

However, this application appears to be a particularly favourable case and it would be of 
advantage to have a larger value of z\> (i.e. larger L) in the matching step to suppress higher- 
order terms. Our analysis suggests that this is indeed possible, since we may determine 
QCD observables rather precisely in, say, the entire range z > 5 by combining QCD results 
at finite z with the static limit. Choosing, for instance, L ps 2Lq in the matching step, one 
needs the QCD observables for z « 10. The errors of the quadratic fit functions evaluated at 
z = z\y = 10 are typically 30%-50% smaller than the errors of the neighbouring points seen 
in figures 6-10. Proceeding in this way, namely taking the fit functions as representations 
of the QCD observables in finite volume (of course within their errors), one may obtain 
the HQET parameters as functions of the quark mass and infer predictions for all quark 
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masses larger than the minimal one considered. With the entire matching done non- 
perturbatively, the final HQET results will differ from QCD by 0(l/m 2 ) if 1/m-terms are 
properly included. No perturbative errors as in Cx enter in this programme. 

We finally remind the reader that in the quantities discussed in this paper, as well 
as in the matching step just described, we are dealing with matrix elements of a mixture 
of energy eigenstates where states with energy of 0(1/L) contribute. Hence the 1/m- 
expansion for the large-volume B-physics matrix elements may be expected to behave even 
better and should be well under control, once this is the case for the matching step. 
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A. Results at finite lattice spacing 



For any details on the lattice simulations and the subsequent analysis of the numerical 
data, the reader is referred to [35] and references therein. In table 2 we therefore directly 
list the results on the observables studied in this work at finite values of the lattice spacing 
as well as in the continuum limit. As already mentioned in section 4, the latter have been 
extracted by linear extrapolations (a/L) 2 — > of the 0(a) improved lattice data following 
the same procedure as adopted in ref. [35]. 



L/a 


YpqfL Ml 


Yv(L Ml 

± V V 1 ) 


R(L Ml 


LAr(L Ml 


E(L Ml 






z 


on 

— o.U 






Iz 


— 1.35(1) 


— 1.581( 7 ) 


U.oby^o J 


— U.oUy(lUj 


n O A K 1 A \ 

0.245(4) 


lb 


— 1.357(7) 


— 1.576(7) 


0.874(b) 


— 0.51z(10) 


0.z5b(5) 


on 
20 


— 1.35/(6J 


— 1.5 /5( / J 


0.8 (0(0 J 


— 0.521(11) 


0.259(6 J 


24 


— 1 357(61 


— 1.576(7) 


874(51 


—0 522(131 


259(81 


32 


-1.360(7) 


-1.577(7) 


0.875(5) 


-0.533(14) 


0.246(9) 


CL 


-1.359(8) 


-1.576(8) 


0.575(e) 


-0.537(15) 


0.250(9) 






z 


= 3.8 






12 


-1.39(1) 


-1.591(7) 


0.884(8) 


-0.419(13) 


0.226(4) 


16 


-1.392(7) 


-1.584(7) 


0.890(6) 


-0.423(12) 


0.240(5) 


20 


-1.391(6) 


-1.582(7) 


0.892(5) 


-0.430(12) 


0.240(6) 


24 


-1.391(6) 


-1.582(7) 


0.891(5) 


-0.431(13) 


0.241(8) 


32 


-1.394(7) 


-1.583(7) 


0.892(5) 


-0.441(15) 


0.228(9) 


CL 


-1.394(8) 


-1.581(8) 


0.893(6) 


-0.444(16) 


0.232(9) 






z - 


= 5.15 






12 


-1.44(1) 


-1.609(7) 


0.904(8) 


-0.315(17) 


0.200(4) 


16 


-1.438(7) 


-1.596(7) 


0.910(b) 


-0.320(15) 


0.215(5) 


20 


-1.437(6) 


-1.592(7) 


0.913(5) 


-0.328(15) 


0.216(6) 


24 


-1.437(6) 


-1.591(7) 


0.913(5) 


-0.327(16) 


0.218(8) 


32 


-1.436(7) 


-1.589(7) 


0.914(5) 


-0.341(17) 


0.210(9) 


CL 


-1.435(8) 


-1.586(8) 


0.915(6) 


-0.343(19) 


0.214(9) 



Table 2: Lattice results on the quantities of this work for different values of the dimensionless 
renormalization group invariant heavy quark mass, z — ML, with L = Lq w 0.2 fm. The full sets 
of simulation parameters can be inferred from tables 1 and 4 of ref. [35], where also the results for 
ir av were published. The quoted errors cover the statistical as well as the systematic uncertainties, 
including those originating from the fact that in the numerical simulations z can only be fixed within 
some finite precision. Continuum limits are displayed in italics. 
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L/a 


y PS (L,M) 


Yy(L,M) 


R(L,M) 


LA r (L,M) 


~(L,M) 






z = 


6.0 






12 


-1.46(1) 


-1.620(7) 


0.913(8) 


-0.270(20) 


0.186(4) 


16 


-1.462(7) 


-1.604(7) 


0.920(6) 


-0.276(17) 


0.203(5) 


20 


-1.460(6) 


-1.597(7) 


0.923(5) 


-0.283(17) 


0.205(6) 


24 


-1.460(6) 


-1.595(7) 


0.924(5) 


-0.282(17) 


0.207(8) 


32 


-1.458(7) 


-1.593(7) 


0.924(5) 


-0.296(19) 


0.198(9) 


CL 


-1.457(8) 


-1.589(8) 


0.926(6) 


-0.297(21) 


0.203(9) 






z = 


6.6 






12 


-1.48(1) 


-1.629(7) 


0.918(8) 


-0.243(21) 


0.178(4) 


16 


-1.477(8) 


-1.609(7) 


0.926(6) 


-0.251(19) 


0.196(5) 


20 


-1.474(6) 


-1.601(7) 


0.930(5) 


-0.258(19) 


0.198(6) 


24 


-1.474(6) 


-1.599(7) 


0.930(5) 


-0.257(19) 


0.201(8) 


32 


-1.473(7) 


-1.596(7) 


0.931(5) 


-0.270(20) 


0.191(9) 




-1.471(8) 


-1.591(8) 


0.933(6) 


-0.271(23) 


0.197(9) 






z = 


9.0 






12 


-1.56(1) 


-1.678(7) 


0.933(8) 


-0.166(31) 


0.149(4) 


16 


-1.530(8) 


-1.633(7) 


0.943(6) 


-0.182(27) 


0.172(5) 


20 


-1.522(6) 


-1.618(7) 


0.948(5) 


-0.190(25) 


0.176(6) 


24 


-1.520(6) 


-1.612(7) 


0.950(5) 


-0.189(25) 


0.180(8) 


32 


-1.517(7) 


-1.608(7) 


0.951(5) 


-0.200(26) 


0.170(9) 


CL 


-1.515(11) 


-1.600(12) 


0.954(7) 


-0.202(30) 


0.177(9) 






z = 


13.5 






16 


-1.639(8) 


-1.713(7) 


0.961(6) 


-0.106(47) 


0.139(5) 


20 


-1.595(7) 


-1.659(7) 


0.967(5) 


-0.123(39) 


0.149(5) 


24 


-1.583(7) 


-1.640(7) 


0.970(6) 


-0.126(37) 


0.157(8) 


32 


-1.578(7) 


-1.631(7) 


0.973(6) 


-0.132(39) 


0.142(9) 


CL 


-1.571(17) 


-1.619(18) 


0.976(9) 


-0.137(66) 


0.145(14) 



Table 2: (continued). 
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B. Perturbative conversion factors between QCD and HQET 



In this appendix we provide details on the numerical evaluation of the perturbative approx- 
imations of the conversion functions Cx{M/A^) (X = PS, V, PS/V, spin) that translate 
the matrix elements and energies obtained in the effective theory to those in quenched 
QCD at finite values of the heavy quark mass. Also the case of the conversion factor 
C mass (M/A^), which relates the heavy quark's pole mass mg to the renormalization 
group invariant quark mass M, will be addressed. 

Let us begin with the definition of the conversion functions for matrix elements of the 
heavy- light currents. Parametrized with the MS mass to*, implicitly defined through 



m Ms( m *) = m * 



(B.l) 



we write them for X = PS and V as 



Cx(m,)= [2M 2 (TO,)] 7( ' /(26o) exp 



f g(m*) 



I 



dg 



7o Xl 



7 X (g) 



(B.2) 



Here, [3(g) = — g 3 &o — g 5 bi + ... is the four-loop anomalous dimension [43] of the renormal- 
ized coupling g(n) in the MS scheme with the leading- and next-to-leading-order coefficients 
bo = ll/(47r) 2 and 61 = 102/(47r) 4 . In eq. (B.2) we have introduced the anomalous dimen- 
sions in the matching scheme 



7 X G?) = V{7o X + 7 1 Y/+7 2 V + ...}• 
At one-loop order they are the universal ones [60, 61], 



(B.3) 



7o PS 



7o V 



1 

"4^2' 



(B.4) 



and at higher order they are related to the anomalous dimensions 7 X > MS (g) of the corre- 
sponding effective theory operator in the MS scheme via 



if 



X,MS . , X 

7 X + 26 cf 



7 2 X = 7 2 X,MS + 46o(c X + ifk) + 2b lC f - 26 [c^ 



Xi2 



(B.5) 
(B.6) 



For X = PS and V, the MS two-loop anomalous dimensions are known from [62, 63, 64] 
and the three-loop ones from ref. [9]. 8 The coefficients 



2 1 

3 4^2 ' 



„ps 



4 1 v 

3 4^' ° 2 



-1 



-4.2 
-11.5 



(4vr 2 ) 2 ' 
1 



(B.7) 
(B.8) 



8 Note that, since in HQET the anomalous dimension of the quark bilinears does not depend on their 
Dirac structure, 7^'^ = 7^'^ holds for all n. 
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originate from the matching of the effective theory operators renormalized in the MS scheme 
to the physical ones in QCD [1, 65, 58], while the term proportional to 

k = ~h (B9 » 

is due to a reparametrization: the matching was originally done at the matching scale 
expressed in terms of the heavy quark's pole mass, toq. Using the perturbative expansion 
for the ratio fn^/niQ [66], 

(B.10) 



TO* 



= l + kf{m*) + ... , 



m Q 

the pole mass can be replaced by to*. 

Next, we turn to the chromomagnetic operator i^^a ■ B?/V If, for the moment, we 
follow the common practice that in the HQET expansion its matrix element is understood 
to multiply the inverse pole mass, 1/toq, the associated conversion function C mag would 
be given by eq. (B.2) with an expansion (B.3) for X = mag and the universal one-loop 
coefficient [67, 68] 

With the two-loop anomalous dimension ,y™ a s> MS j n the effective theory given in [51, 52] 
one finds 

7 mag = 7 mag,MS + ^mag ^ ^ = lj! JL , (B.12) 

where c™ ag was determined in [67]. In view of eq. (3.19), however, we are rather interested 
in a function C spm , which multiplies 1/M. In other words, the conversion function C sp i n 
must also include the factors to*/to.q and M/to* in order to cancel the conventional factor 
1/toq in the HQET expansion in favour of 1/M. Using the relation 



M 



— = [2b g 2 (m*)] 

TO* L J 



-do/(26o) 



exp < 



/ 



dg 



rjg) _ _^o_ 
0(g) b g 



(B.13) 



where r(g) = —g 2 do — g A d\ + . . . denotes the quark mass anomalous dimension in the MS 
scheme in QCD (known up to four loops [45, 46] with leading coefficient do = 8/(4-7r) 2 ), 
and the ratio (B.10), we then obtain eqs. (B.2) and (B.3) for X = spin and 

6 



spin 

7o 



mag 

7o 



d 



(4vr)< 



- d , 



spin mag , n , , mag, MS , . t / mag . ; \ 

7i = 7i -di + 2b k = j 1 6 ' -di + 26 (c 1 s + k) , 



(B.14) 
(B.15) 



with ^/™ ag from eq. (B.12). 

For the special case of the ratio of pseudoscalar and vector current matrix elements, 
X = PS/V, all but the contributions from the matching cancel and one gets 



Cp S /v(to*) = exp < 



g(m*) 



dg 



7 Fb G?)-7 V (ff) 



(B.16) 
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To parametrize the mass dependence of energy observables (such as T av in eq. (3.12) 
of section 3.2), we also define 



r t- \ - m ® m * 
m* M 



(B.17) 



where in this case the highest available perturbative precision is achieved by taking the 
four-loop r-function [45, 46] in rn*/M together with the three- loop expression for rag/m» 
[47]. 

Finally, changing the argument of the various Cx to the renormalization group invari- 
ant ratio M/Aj^g via (B.13), we straightforwardly arrive at the conversion functions 

C x (M/A m ) = 6 X (m*) with X = PS, V, PS/V, mass, spin . (B.18) 

We evaluate all occurring integrals in the above expressions exactly (sometimes numeri- 
cally), truncating anomalous dimensions and the /3-function as specified. Also eq. (B.l) is 
solved numerically. For practical purposes, such as repeated use in the fits of the heavy 
quark mass dependence of our QCD observables, a parametrization of all conversion func- 
tions in terms of the variable 

x = , - r- (B.19) 

In (M/A m ) 

was determined from a numerical evaluation. The functions decompose into a prefactor en- 
coding the leading asymptotics as x — ► 0, multiplied by a polynomial in x, which guarantees 
at least 0.2% precision for x < 0.6: 

( ~,PS/(9hrA r -, ^ n„r- n„,n 91 ^ , PS 



C PS (x) 



x 7 7(2&o) | i _ 0.065 x + 0.048 a; 2 } 



2-loop 7 



k x^o 7(2&o) | i _ 0.068 x - 0.087 x 2 + 0.079 x 3 } 3-loop 7 



, (B.20) 



PS 



C v (x) 



x^ /( 26() ) { 1 - 0.180 x + 0.099 x 2 } 



2-loop 7 



v 



k x^ /( 26 °) { 1 - 0.196 x - 0.222 x 2 + 0.193 x 3 } 3-loop 7 V 



, (B.21) 



' 1 + 0.117 x - 0.043 x 2 



2-loop 7 PS ' V 



1 + 0.124 a; + 0.187 x 2 - 0.102 a; 3 3-loop 7 PS ' V 



(B.22) 



' a;<V(2M I 1 + o.247x + 0.236 a; 2 } 



C mass (x) — < 



2-loop 



x do/(2b ) | 1 + .i79a; + 0.694a; 2 + 0.065 x 3 } 3-loop 



, (B.23) 



fn-t 



x 7o pm /(2feo) { 1 + 0.O66 x } 1-loop 7 spin 

x 7o pm /(2bo) I 1 + o.087 x - 0.021 x 2 } 2-loop 7 spin 



(B.24) 
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Apart from the function C mass , the pole mass does not appear in any of the above per- 
turbative expressions; they relate observables in the effective theory to those in QCD and 
are parametrized by the RGI mass M, which is unambiguously defined in terms of the 
(short-distance) running quark mass (see eq. (3.6)). Thus, their perturbative expansion 
is expected to be a regular short-distance expansion. In particular, it is not expected to 
suffer from the bad behaviour of the series (B.10). 
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